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. Abstract 
^ \ Self Consistent Quasiparticle Random Phase Approximation (SCQRPA) is considered 

^ ■ in application to the Fermi transitions within the 0(5) model. It is demonstrated that 

^ ! SCQRPA improves on renormalized QRPA (RQRPA), a method that has recently become 
■ - - ' rather popular in this context. The analytical form of the SCQRPA vacuum is used to 
evaluate all the matrix elements. The SCQRPA results show a general trend similar 
to the exact solutions. The necessity to change the single particle basis beyond the 
transition point, and to include the proton-proton and neutron- neutron channels in the 
QRPA operator, in addition to the proton- neutron one, is pointed out. 
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1 Introduction 



Since Vogel and Zirnbauer and Cha have discovered the importance of the particle- 
particle force in the S = 1, T = channel, the quasiparticle random phase approximation 
(QRPA) has been the most frequently used nuclear structure method for evaluating double 
beta (PP) rates both for the two-neutrino decay mode {PP2u) and for the neutrinoless 
mode {f3(3ou)- The general feature of this method is that the resulting nuclear matrix 
elements A42u and Aioi, turn out to be highly sensitive to the particle-particle force 
|jl|, I], ^, |, ^, 1^, |[. Moreover, very close to its physical value the QRPA collapses. One 
thus may suspect that this method yields relatively small values of A42u simply because 
the approximation breaks down. In other words, the smallness of Ai2u in the QRPA could 
be just an artifact of the model. 

In the past few years, several modifications of the QRPA have been proposed to 
improve on the above mentioned behaviour |10|, |Tl], [l^, Yet, the one that has 
received major attention lately is the so-called renormalized QRPA (RQRPA) [|14[ |15|, 
T7|, p!8|, |19|, |20, pl| . The new ingredient that is brought up is the effect of the ground state 



correlations (GSC) in the QRPA equations of motion (EOM) i.e. while in the standard 
QRPA the correlated ground state is approximated by the BCS vacuum, in the RQRPA 
the GSC are partially taken into account. Let us recall that the GSC are fully incorporated 
only in the self-consistent RPA (SCRPA) theory |2^, which simultaneously leads to a 
coupling of the single-particle field to the RPA excitations. The same formalism was also 
generalized to the superconducting systems , and when applied in this context will be 
referred to as the self-consistent QRPA (SCQRPA). 

The RQRPA generally yields better results than the QRPA, in the sense that the 
instability is avoided, and the matrix elements A42u and J^Of are somewhat less sensitive 
to the particle-particle interaction. Yet, the charge-exchange sum rules both for Fermi 
(F) and Gamow- Teller (GT) transitions, are significantly violated within the RQRPA 



It has been shown [|1^ that the origin of this violation is the omission of the 
contribution of the so-called scattering (anomalous) terms to these sum rules. A possible 



where the inclusion of these 



way to solve this problem has been suggested in ref. 
terms in the definition of the excitation operators is proposed. 

To explore the effect of the GSC on the f3f5 matrix elements, several applications of 
the RQRPA have recently been performed in the framework of the exactly solvable 0(5) 
model which allows for a direct comparison between the exact and approximate 

results. A similar exercise has also been done within the SCQRPA ||28|, but by solving 
the corresponding EOM in an approximate way. For that reason, and in contrast with 



experience with the SCQRPA approach |2^, | 
elements for the PP decay were not well reproduced. 
In this paper we resume the study performed in ref. 



the excitations energies and matrix 
28| , but avoiding the approxima- 
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tions done there. With the aim to clarify the physical meaning of the model parameters, 
a brief introduction of the 0(5) model is presented in Section 0. The full SCQRPA equa- 
tions for the 0(5) model, with the simplifying condition of decoupling the charge-exchage 
mode from the conserving ones, are derived in Section ^. These equations are solved 
self-consistently with the generalized BCS equations, which implies the coupling of the 
quasi-particle mean-field to the charge-exchage excitation. The task has been facilitated 
by the circumstance that it was possible to construct explicitly the correlated ground state 
wave function for the schematic 0(5) Hamiltonian. In Section ^ the RRPA is formulated 
as the limit of the SCQRPA when the ground state two-body density is approximated 
by products of one-body densities. Numerical results for the Fermi transitions and the 
corresponding (3(3 decay matrix elements, obtained within different approximations, are 
compared with each other and confronted with exact results in Section ^. Summarizing 
conclusions are given in Section ||. Finally, some details of the calculations are displayed 
in the appendices. 



2 The Model 

In the model, protons and neutrons occupy only one level and the Hamiltonian 



H 



+ 4 



1 



G 



pp\p 
^ ^pp^pp 



G 



p (jv jv +/V jv ) - G (P^ P + P"^ P 

pn V-' 'pn-' 'np ' ^ "np-'- " pnj ^pn \^ pn^ pn ' ^ np np 

is such that it can be written in terms of generators of 0(5) |2 

5Z ^am^bmi 



N. 



ab 



plb = n!, = E(-i)^'-^c:^c, 



am b—m^ 



(2.1) 



ab 



Pi 



ab I ' 



(2.2) 



where ab stands for pp, nn, pn and np. 

To see the physical meaning of the model parameters, we rewrite the model Hamilto- 
nian ( p. ID in the coupled basis 



^ab:JM 



(4 ® 4) 



JM 



(2.3) 



I.e. 



H 



a=p,n ^ t: / ^ a=p,n 



aa-m 



2 



t -^P" J t 

^GpnCpn-ooCpn-00 -\ 7^^^{~^) Cpn;J M^pn-J M (2-4) 
^ JM 



and compare it with the standard one-level shell model Hamiltonian: 

- J-^aa + J 2^ [i + (-ij J , ^aa;JM^aa;JM 

a=p,n ^ JM;a=p,n ZJ -\- L 

, {jpjn, J\ \v\\jpjn; J) + {-lyjjpjn, J\ \v\ Ijnjp] J) n (O 

+ / o T I 1 '^pn;JM'^pn;JM- [^.0) 

JM V^-' + i 

We conclude that the e^^ are related with ea as 

C = + ^, (2.6) 

and that the particle-particle matrix elements for the Hamiltonian are: 

G{aaaa- J) = [1 + (_i)^]Mij^Hb|Z) = -2fiG,5jo, (2.7) 

and 

{jpjn, JMljpjn, J) + (-l)'^(jpjn; J\\v\\jnjp; J) 



G{pnpn; J) = 



V2J+1 



with 2Q = 2j + 1. Therefore like particles interact only in the J = channel (pairing 
term). But, the unlike particles, besides having a pairing term, also interact in J 7^ 
channels. In particular, for odd J, the Hamiltonian ( |2.1| ) simulates the T = pairing. 
Moreover, it is charge dependent, and the limit of charge independence is achieved when 



dp 6yx, Gp Gji GpYi and Fpn 0. 



For further use it is convenient to define the particle-hole matrix elements Q 

F{aaaa; J) = -(-l)-^Ga (2.9) 



^The particle-hole matrix elements are defined as |33 : 

F{adcb;J)^-J2i'2J' + ^)[l \ J } 



G{ahcd] J'). 
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F{pnpn; J) = (qF^Jjo - (-1)"^) , (2.10) 

F{ppnn; J) = [f^^ + {-IYG^^) . (2.11) 

It might be worth noting that to make the Hamihonian ( |2.1D charge conserving, even 
when Fpn 7^ 0, one has to do the replacement: 

H^H+ ^NppNpp + ^NnnNnn. (2.12) 

Now 

G{aaaa-J) = -2G',n5jo + y (1 + (-1)^) 

F{aaaa;J) = -(-1)-'G„ + y (2fi5jo - 1) (2.13) 

and the charge independence is fulfilled for = e„, Gp = = Gp„ and Fp = Fn = Fpn. 
In this case the following relations are valid: 

G{aaaa]0) = 2G{pnpn]0), 

F{aaaa;0) = F{pnpn;0) + F{ppnn;0). (2.14) 



3 The SCQRPA Equations 
3.1 General Formulation 

The "nuclei" to be considered here are usually superfiuid both for protons and for neu- 
trons. We therefore make first a Bogoljubov transformation to quasiparticle operators 
at 

"'ami 

cim = Uaaim + Va{-iy^'^aa,-m, a = p,n. (3.1.1) 

The SCQRPA states are constructed by the action of the excitation operators 

rt= X^^Al-Y:Ar, //= 1,2,3 (3.1.2) 

T=pp,nn,pn 

on the SCQRPA ground state |0), 

k)=rt|0), (3.1.3) 

^ In recent works [ p3| , it was proposed to also include anomalous terms a^a in the excitation 
operator. Though this potentially is of great importance we do not consider this extension here (see also 
the discussion at the end of this paper). 
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with the vacuum condition 
The operators 



r,|o) = 0. 



ab 



(3.1.4) 
(3.1.5) 



I2n{l + dab) 

which appear in eq. (|3.1.2| ) create pp, nn and pn quasiparticle pairs with J = 0: 

Al = E(-l)''""«L<-m, = P^n. (3.1.6) 



and therefore the ansatz ( p.l.2|) can be considered as a Bogoljubov transformation between 
fermion pair operators. 

The normahzation condition 



{^W) = (o| 



|o) = s, 



vv' 1 



gives 



where 



J2 - y^^:*d;,^y^^ 



D,>, = {0\[Ar',Al\ |0). 
The SCQRPA equations have then the form: 

[DO 



A B 

B* A* 



X 
Y 



-D* 



X 
Y 



(3.1.7) 
(3.1.8) 

(3.1.9) 
(3.1.10) 



where 



A,,. = (o| |o) = ^;v, 

Brr' = -(0| [Ar,n,Ar'] |0) = B^,^, 

are the symmetrised double commutators defined by Rowe ||3^ with 



Ti — H — XpNpp — XnNnn- 



(3.1.11) 



(3.1.12) 



The SCRPA equations are still incomplete because the mean field parameters Ua, Va, 
and Xa are so far undetermined. This general problem of SCRPA has been solved very 



satisfactorily in where it has been shown that the mean field parameters (the single 
particle basis in a non-superfluid system and/or the amplitudes Ua, Va in a superfluid 
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one) can be determined in minimizing the SC(Q)RPA ground state energy with respect 
to these parameters. This minimization turns out to be equivalent to p4| : 



(0|[7^,rt]|0) = 0; z/=l,2,3, 



(3.1.13) 



which with ( p.l.4|) is a very natural equation from the EOM point of view. 
The eq. ( |3.1.13| ) is equivalent to 

(0|[7^,4]|0)=0, 



(3.1.14) 



which, upon replacing |0) by the mean field ground state, turns out to be identical to the 
usual BCS equation. Q 

The chemical potentials \n,p are determined from the usual particle number condition, 
i.e. it is imposed that the SCQRPA ground state |0) has, on average, the correct number 
of protons (Z = Np) and neutrons (A^ = N„) 



(0|Ar,jO) = N„ 



(3.1.15) 



The set of equations ( p.l.4|) , ( |3.1.1(]| ), (|3.1.11|) and (|3.1.15| ) have to be solved simulta- 
neously, together with the conditions (|3.1.14|) . This generally makes the problem to be 
very demanding computationally. Indeed the problem has only been tackled so far within 
very simple schematic models [|^. Here we are going to work out the SCQRPA equa- 
tions of the 0(5) model, under some simplifying conditions. That is, we assume that the 
charge-exchange mode. 



rl 



pn 



X, 



pn 



pn 



A. 



■pn 



/2n '"^/2<a' 

decouples from the modes involving like particles pairs. 



Pt 



E 

a=p,n 



2vn 



A 



2vn 



(3.1.16) 



(3.1.17) 



Since both the Hamiltonian and the Bogoljubov transformation are charge conserving, the 
decoupling happens automatically for a SCQRPA vacuum that is a superposition of states 
with even number of protons and neutrons. This is achieved by taking into account only 
the correlations due to the pn pairs. In other words, the SCQRPA vacuum is determined 
only by the equation 



PpnlO) = 0. 



(3.1.18) 



■^These generalized mean field equations also serve to make the double commutator s (3.1.1l| ) symmet- 
ric. In the 0(5) model the symmetry of the double commutators is verified (see eqs. (3.5.9)). 
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Within this framework, and since 



pn 



(3.1.19) 



is trivially satisfied, the eqs. ( |3.1.15| ) for the mean field reduce to 

{o\{n,AlM = {o\[n,AU\o) = o, 



(3.1.20) 



and the number condition (|3.1.14| ) has to be obeyed. 

In summary, since in this paper we are interested in investigating the Fermi transitions 
in this schematic model, our equations are: (|3.1.16|) for the excitation operator, (|3.1.18|) 
for the SCQRPA vacuum, and ( |3.1.14| ) and ( |3.1.19| ) for the mean field. 



3.2 Hamiltonian in the quasiparticle basis 



The first step in solving the above equations is to rewrite the Hamiltonian ( |3.1.12 ) in the 
quasiparticle basis. Using the Bogoljubov transformation, 7i can be put in normal order 
with respect to the quasiparticle vacuum. Five type of terms occur. 



(3.2.1) 



and they can be expressed in terms of the generators of an 0(5) algebra J^^^, Aah = (-^L) ^ 
and 

Mab = '^al^abm, (3.2.2) 



analogous to those given by ( |2.2| ), but now in the quasiparticle basis. The first term 



a=p,n 



- Q{Gpn + Fpn)VpVl, 

is the BCS ground state energy, and the second one 
is the one-quasiparticle Hamiltonian, where 

Ea={e^^-\a){^il-vl)+2nvlulGa, 

are the quasiparticle energies, in the BCS limit. Here 



(3.2.3) 



(3.2.4) 



(3.2.5) 



(3.2.6) 
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are the single particle energies, corrected by the self-energy terms |]T7 

1 



-Gpvl 



The remaining quantities are 



n 



40 



a 

~l~ 'o^i.^pn ~l~ FpnjUnVnUpVpi^^^yi^^^pj^ -\- •/A.pn-^pn) ■, 



and 



n 



31 



n 



22 



^ a 



a a 



+ -^[FpniUpU^ + VpVj - GpniUpV^ + U^V^)] 
X ij^pn-^np ~l~ •^np^pn ■^nn -^pp)- 



3.3 The SCQRPA vacuum 

The correlated vacuum is defined by eq. ( ^.1.181 ), which is equivalent to 

Y 



Mp„-2;4J|0)=0 
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From the vacuum condition it is very difficult to find an explicit expression for the SC- 
QRPA ground state |0) in the general case. Yet, because of the simplicity of the model 
used here, the eq. ( p.3.1|) can be solved exactly, and one obtains 

|0) = ^/iVoEai;^' {Al^f \BCS) = ^oY.aiz'\2l), (3.3.2) 

1=0 1=0 

with \BCS) being the quasiparticle vacuum and 

N,= {j2a',z'\2l\2l)^ , (3.3.3) 

, ^„ Q\ (2Q-2l))\ ,,,,, k\(2Q)\ 

Thus, the SCQRPA vacuum is a superposition of states with equal even number of neu- 
tron and proton quasiparticles, and as such a superposition of states with even number 
of protons and neutrons. To solve the SCQRPA equations we have to evaluate the expec- 
tation values of the generators A^^,, Aab and J\fab and their bihnear combinations as well. 
The results are given in the Appendix A. 



3.4 Generalized BCS equations 



The only undetermined quantities are the u^s and v^s of the BCS state and the chemical 
potentials. They are fixed by the generahzed BCS equations (|3.1.19|) and the number 
conditions ( p.l.l4| ). The first ones yield 

{o\[n'' + n'\Al]\o) = 0, 

{o\[n'' + n'\ALm = 0, (3.4.1) 

and from the second ones we obtain 

Na = 2Qvl + {ul-vl){0\Ka\0), a = p,n. (3.4.2) 

When the SCQRPA vacuum |0) is replaced by the quasiparticle vacuum \BCS), the eqs. 
( ^.4.11 ) and ( |3.4.2D reduce to the usual BCS equations 



2(6^^) - A,) - - 
2(eW - A„) - GM< - 



UpVp 



0, 

UnVn = 0, 



(3.4.3) 



{BCS\[n''^,Al]\BCS) = 

{BCS\[n^'^,AlJ\BCS) = -zur -Ar,)-u„^i[uz-v: 

This generalized particle number equation has for the first time been incorporated into an extended 
RPA calculation in refs. Q though this had not explicitly been mentioned in the latter of the two 
references. 
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with 



2fi' 



p, n. 



The SCQRPA correlations lead to the following generalized BCS equations: 



0; 



2 I 2 



(3.4.4) 
(3.4.5) 



with the standard solution for u's and v's 

n 1 



k= 



(3.4.6) 



where 



\p){om-Afpp)\o) 

^ {Fpn + Gpn) {vl - ul) (0| {Al^Apn ' MnpMpn) |0), 



(3.4.7) 



and 



A. 



2VL 
1 



GpUpVpi<d\ [2(^] -Tv; 



A 



|0) 



— {Fpn + Gpn) UnVn{0\Al^AUO) 



2n 



(3.4.8) 



are, respectively, the renormalized single particle energy and gap parameter for protons 
(and similarly for neutrons). One remarks that formally the generalized BCS equations 
have the same structure as the ordinary ones. However proton and neutron equations are 
coupled, the pairing interaction is renormalized by the two body densities (screening) and 
also the single particle energies are coupled back to the correlations. 

From the number of particle condition (|3.4.2| ) it follows immediately that the BCS 
amplitudes ( |3.4.6| ) can be expressed as: 



1 
2 



n-N„ 



1 



n-{o\Ka\o), 



It should be stressed that these expressions, which generalize ( p.4.4|) are independent of 
( p.4.6|) and are equivalent to the eq. ( p.4.2|) , which fixes the number of particles. They 
allow us to express directly the quantities entering the BCS eqs. ( p.4.5|) in terms of the 
one-body densities and to solve them. 
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3.5 The SCQRPA equations 

As observed previously, the vacuum state |0) is a superposition of states with even number 
of protons and neutrons. Therefore, since the Hamiltonian is charge conserving, the 
SCQRPA eqs. (|3.1.1CI| ) spht into two equations, one for the charge-conserving excitations 
and the other for the charge-exchange excitations. Being interested only in investigating 
the Fermi transitions, we will consider only the last ones. 



n 



pn 



pn ' -ni 



/2n ""^J2n^ 
\pn) = rt |0). 



with 

The normalization condition 
gives 



(pnM = (o|[r,„,rt ]|o) 



Next, we introduce the normalized pair creation operators 

-1/2 -^pn 



A^ 



pn 



and the normalized amplitudes 



pn 



^pn ^pnj 



J'pn -^nn ^ I 



pn ^pni 



with 



\Apni A'L 



pnl 



to write the SCQRPA equations as 



|0) 



{Q\Mpp + Mnn\Q) 



2VL 



A B 

B* A* 



■pn 



UJ 



pn 



1 
-1 



pn 



pn 



Here u is the SCQRPA excitation energy and 

A = {0\[A,nAn,Alj]\0) = {Q\[[Apn,n],AlJ 



|0), 



B 



(3.5.1) 

(3.5.2) 

(3.5.3) 
(3.5.4) 

(3.5.5) 

(3.5.6) 
(3.5.7) 

(3.5.8) 



(3.5.9) 
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Using the expression for H in the quasiparticle basis, one can evaluate the double com- 
mutators in eq. ( p.5.7|) , and obtain 



A 
B 



^11 ^ p2 ^ 



40 



^40 ^ ^22 



(3.5.10) 



with 



A'^ = {o\[A,nAn'\Al]\\o) 

B'^ = -{0\[Apn,[>i'',An]\\0). 
The explicit expressions for A^^ and B^^ are: 



A^^ = Ep + En 



(3.5.11) 



(3.5.12) 



with Ea given by eq. (|3.2.5|) , and 

A'' : 



(0|4„^p„|0) 



+ [Gp{ul + vD + GnK + O + Fpnfi - Gpna 
- A{Gpulvl + GnU^V, 



(0| (2fi- AT) |0) 

{Q\AlpApp-2UpnUnp\^) 



(3.5.13) 



(0| (2f]-Ar) |0) 



2 2 , r. .,2..2^(0|(2^^-A^)A^|0) 
IpVp 



(0|(2f]-Ar)|0) ' 



40 



iGpn ~\~ Epn)UpVpUnVn 



{^\ApnApn ~^ ApnApn\0) 

^)K2?r~A0|0) ' 



B 



40 



'{Gpn -\- Fpn)UpVpUn'Vn I 1 



{2n-Aff-2AlApn\0) 



(0| (2fi-Ar) |0) 



(3.5.14) 
(3.5.15) 



B 



22 



[Gp{u +v ) + Gniun + t^„) + FpnP - Gpna 



{0\ApnApn\0) 

(o| {2n-M) |o) 

(0|^ppA„|0) 



(0|(2r]-A/')|0)' 



where the short notations, 



TV = Mn.n + Mr 



pp I J V } 



(3.5.16) 
(3.5.17) 
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and 



2 2 I 2 2 



(3 = uX + vyp, (3.5.18) 
have been used. 

It easily can be seen that, when the ground state correlations in the SCQRPA equations 
are neglected, i.e. when the SCQRPA vacuum |0) is replaced by the quasiparticle vacuum 
\BCS), the usual QRPA equations 



A = —{Gp + Gn)—/3G{pnpn;0)+aF{pnpn:,0), 
B = [F{pnpn; 0) — G{pnpn; 0)]unVnUpVp, 



(3.5.19) 



are recovered. 



4 Renormalized QRPA equations (RQRPA) 

As pointed out before, to solve the SCQRPA (SCRPA) equations exactly is very demand- 
ing from a computational point of view, and this has not been done so far in realistic cases. 
Several authors ||T^ |T^, |T^, |T^, [T^, ^ |HT|, ^ have suggested different approximations 
to derive simpler equations to be amenable for the numerical calculations. In essence, the 
simplifying hypothesis consists in factorizing the two-body correlations in the SCQRPA 
(SCRPA) equations of motion. For a normal system this implies to approximate the 
ground state two-body density by products of one-body densities. 



(0|cHc,c,|0) = (0|clc,|0)(0|4c,|0) - (0|clQ|0)(0|4ce|0). 
The parallel approximation in a superfluid system is: 



(4.1) 



{0\alaladac\0) 




(0 


alac 


0)(0 




0)- 


(0|atGd|0)(0|4a,|0) 




+ 


(0 


4«i 


0){0 




0), 




(0|al44«c|0) 




(0 


4«c 


0)(0 




0)- 


(0|44|0)(0|4a,|0) 




+ 


(0 


44 


0)(0 




0), 




(01444410) 




(0 


44 


0)(0 


aWd 


0)- 


(o|44|o)(o|44|o) 




+ 


(0 


t t 


0)(0 




0), 





(4.2) 

and the corresponding hermitian conjugates. 

The SCQRPA ground state does not have a fixed number of quasiparticles. Thus its 
quasiparticle densities have both terms that conserve the number of quasiparticles and 
terms that violate it. In field theory language, one would say that in eq. ( p.5.1|) we have 
approximated the four-point function by a product of two-point functions. Note that 
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even with this approximation, to calculate the one-body densities we have to evaluate the 
ground state |0) explicitly, and this is very difficult in practice. Therefore, the one-body 
density (two-point function) is approximated so that it is not required to determine the 
excitations and the ground state simultaneously. This class of approximations are known 
as RQRPA (RRPA). They are valid as long as the effect of the two-body correlations can 
be neglected, and their main advantage, in comparison with SCQRPA (SCRPA), is the 
numerical simplicity. 

In this paper we derive the RQRPA equations for the 0(5) model following the scheme 
described above. The RQRPA limits for the expectation values are shown in the Appendix 
B. Within this limit, the generalized BCS equations (|3.4.5| ) reduce to 



ASM) _ ^ 

Cp /\p 



Gr. 



- -{Fpn + G. 



pn ) 



2n 



0, 



(4.3) 



and an analogous equation for neutrons. 

In the BCS limit the self-energy terms for protons are given by (|3.2.71 ), i.e. by 
F{pppp;0) = —Gp times Np/2f2, plus F{ppnn;0) = —{Fpn + Gp„)/2 times Nn/2Q. In 
the RQRPA one has to consider the average number of neutrons and protons in the SC- 
QRPA ground state |0) (see eq. ( p.4.2|) ). Thus, the self-energy terms in the RQRPA have 
the same interpretation as in the BCS case. 

Using the number condition (|3.4.2| ), the eq. (|4.3| ) can be rewritten as 



4'"' - A, - Gp- - (F, 



pn 



Gp^Dpp ^ _ ^ 
"4^2 2 ^ ^ 



0, 



(4.4) 



with 



1 - 



{0\Ka\0) 



(4.5) 



Thus in the RQRPA limit the generalized BCS equations are equal to the usual ones, 
except for the renormalization of the gap parameter by the factor Dpp. 

On the other hand, in this limit A gets contributions only from the 11 and 22 terms, 
i.e. A = A^^ + y4^^, and B only from the 40 term, 2. e. B = B^^. One obtains 



^ = -^(GpDpp + GnDnn) - l3G{pnpn; 0)Dpn + aF{pnpn; 0)Dp„, 
B = [Fijpnpn; 0) — G{pnpn; 0)]DpnUnVnUpV. 



pi 



(4.6) 
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where a and f3 are given by ( |3.5.18| ) and Dpn by ( p.5.7| ). Thus, we see that, compared 



with the QRPA equations (|3.5.16|) , the matrix elements F{pnpn;0) and G{pnpn;0) are 
now renormahzed by Dpn, as aheady was known previously |jT^. Yet, our careful renor- 
malization procedure shows that the BCS quasiparticle energy 

BCS ^ ^'^g ^ G{aaaa; 0) 

2 4 ' ' ) 

is renormahzed as 

= ^f^^^aa- (4.8) 

This renormalization is a cooperative effect coming both from A^^ and A^^. In particular 
the former is a consequence of the coupling of the mean field to excitations, embodied in 
the generalized BCS equations 



5 Numerical results and discussion 

An exhaustive numerical study of the 0(5) model has been performed for different values 
of N, Z and Q, as well as for a large set of parameters. For the discussion we have selected 
the following three samples: 



(A) 
(B) 
(C) 



n = 5, N = 8, Z = 2, Gp = Gn = OA MeV, F„„ = 0.4 MeV 



p "^71 xvj-v^ Y , J. pn 

n = 10, = 14, Z = 6, Gp = Gn = 0.3 MeV, Fp„. = 0.3 MeV, 
n = 25, N = 30, Z = 20, Gp = Gn = 0.2 MeV, Fp„ = 0.2 MeV. 

The strength Gpn, or more precisely the ratio, 

Gr, 



' pn 



Gp 



5 



To generalize the result (18) to the many level case, where 



.BCS _ _ 

AUaVa ' 



— {2ja + V(2jb + l)^'\ibVbG{aabb; 0), 

there are at least two possibilities: 

Dg 
^UaVa ' 



= -7^(2ja + 5](2j, + l)^^'u,v,G{aabb; 0), 



and 

4:UaVa 



e"^' = -^^{2ja + V(2jb + lY/^UbVi.G{aabb- 0)Dab 

„7!„ ^ ^ 



(5.1) 



The last one has been obtained in ref. [ p7[ from the self-consistency between the residual interaction and 
the mean field. See also ref. iBOl. 
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is kept as free parameter. The values of the single-particle energies are of no relevance in 
the context of the present work and we have simple chosen ep = e„ = 0. We believe that 
the last case is as "realistic" as it can be within a schematic one-level model. Namely, 
we mimic the calculation of a nucleus with A = 50, within a rather large configuration 



space, with the Kisslinger and Sorensen |Q estimate for the pairing coupling constant 



(G = 25/ A MeV), and the strength Fpn of the order of magnitude of Gp {Fp^/Gp = 1) 
5.1 Solutions of the RPA equation 



By solving the equation ( |3.5.8| ) we get the frequencies u and the amplitudes X and 3^. 
From the last ones we can evaluate the (3^ strengths for the transitions from the even-even 
nucleus (A^, Z) to the odd-odd nuclei {N ± 1, Z ^ 1), i.e. 

s^ = j:m\T^m'; (5.2) 



where T+ = Nnp and T_ = Npn. In an exact calculation the summation in eq. (|5.2|) goes 
over all states in the odd-odd nuclei. But, within a QRPA-like calculation done here there 
is only one intermediate state \0,„t), and 



(0„,i|T_|0) = J2ilDp„{XupVn + yvpUn) 



(0,„,|T+|0) = j2nDpn{yUpVrr + XVpUn). (5.3) 



It has been pointed out in ref. pO] that the sum rule 



S' - = N - Z, (5.4) 



is violated in the RQRPA. However, it was shown |^5[ [T^, ^ that this does not happens 



when both the correct particle number condition ( |3.4.2| ) is used and only J = ground 
state pn correlations are included. Under these circumstances (0|A^|0) equals (0|A/'„|0) 
and the above sum rule is fulfilled in RQRPA as well as in SCQRPA. 

In figures 1 and 2 are plotted the numerical results for uj and S~ , respectively. We do 
not show the results for 5*"^ as they follow from 5*" and the sum rule ( |5.4| ). As expected, 
for UJ the RQRPA results always fall in between those of the QRPA and SCQRPA. The 
QRPA collapses at a value of s = Scrit that decreases when Q is increased (scrit = 2.10, 1.41 
and 1.10, for the cases (A), (B) and (C), respectively). On the contrary, in RQRPA and 
SCQRPA there occurs no collapse and u approaches zero asymptotically when s oo. 
From the figures we see that for s < Scrit, all three approximations agree, but near the 
transition point and beyond they differ. This behaviour is more evident in case (A). The 
exact results for the transition strengths S~ are also shown in figure 2. As evidenced from 
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this figure, for s ^ Scrit tlie approximate S"^ strengtlis coincide so well with the exact ones, 
as well as among themselves, that is hard to distinguish them visually. When s approaches 
Scrit the QRPA values of increase very fast and go to infinity. The SCQRPA results also 
reproduce nicely the exact results in the neighborhood of the critical value of s in all three 
cases. Moreover, in the case A, SCQRPA yields approximately correct strengths even 
for values of s very far from the physical region for this parameter (s = l).The variations 
of z, {Af) and in the case C, as a function of s, are illustrated in figure 3. We see that: 
i) the approximate calculations of the number of quasiparticles in the ground state of the 
even-even nuclei only differ from one another for s > Scrit, and ii) they undergo a sudden 
change for s ~ Scrit, where the QRPA value goes to infinity. The variation of v"^ with s in 
SCRPA and RQRPA is a direct consequence of the particle number condition (|3.4.2|) and 
is more pronounced in the latter case. 

5.2 Excitation energies of the odd-odd nuclei 

While the RPA results for transition strengths S"^ can be directly compared with the cor- 
responding exact values, as well as with the experimental data in a realistic case, this does 
not happens with the pn-QRPA energies u. In fact, one important question is to establish 
the relationship between u and the excitation energies of the intermediate odd-odd 
nuclei {N ± 1, Z ^ 1) from the even-even nucleus {N,Z). This question is particularly 
relevant for the decay. Actually, in most of the related theoretical studies done in the 
framework of the QRPA, this problem is circumvented by the use of the experimental 
ground state energies of the initial, intermediate and final nuclei, and only the relative 
spacing of the excitation energies is coming from the roots of the QRPA equation. (See, 
for instance, ref. 



In ref. 0] (see also ref. [0) it was proposed that the relationship between u and 
is given by 

E^ = ujTXp± A„, (5.5) 



in full analogy with the simple BCS case |35 



However, from a more fundamental point of view, in the framework of the SCQRPA, 



it is more natural to define these excitation energies as |^| 



E^^ = u;t mVnp, [XpNp + A„Ar„,rtJ]|0), (5.6) 

which leads to 

E^^=ujTiX' + y')j:>^aiul-vl). (5.7) 

a 

To test the quality of these two prescriptions we should remember the following prop- 
erties of the energies 
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1) for ep = e„ and Fp„ = 0, E~ should pass through zero at s = 1, and 

2) when the energy shift A = — e„ (which can simulate the Coulomb energy displace- 
ment) is introduced, =p A. 

Obviously, in the exact calculation both of these conditions are fulfilled. Additionally, 
from the eq. ( |3.4.5| ) and the fact that u does not change by a global shift of the single 



particle energies, we see that while the relation (|5.5|), being linear in Ap — A„, obeys the 



condition 2), the prescription (|5.7| ) does not. On the other hand, from the results shown 
in figure 4 for the case (B) with Fpn = 0, we also conclude that only the first prescription 
for E~ is consistent with the condition 1). Thus, from now on we will not discuss anymore 
the energies E^ based on the eq. ( |5.7| ). Why this last definition of the excitation energies 
leads in our context to unphysical consequences will be investigated in a future work. 

In figure 5 are compared the approximated energies E~ with the exact ones. All three 
approximations exhibit quite similar behavior for s^Scru- Within this range of values for 
s, they also agree reasonable well with the exact results. Yet, for s > Scrit the agreement 
between the exact and the RQRPA and SCQRPA results is not anymore so good and 
worsens when Q is increased. In any event the SCQRPA is preferable to the RQRPA, 
both close to the critical point and beyond. 



In all previous studies of the 0(5) model p6i P^ , [28[] only the excitations in the 
(iV — 1, Z + 1) systems have been studied. But, the energies E~^ could be as interesting for 
examining the different approximations as are the former ones. Thus, we show them in 
figure 6. The first issue that attracts attention, when comparing with the results presented 
in figure 5, is that the discrepancies among different QRPA theories for E~^, in the cases 
(A) and (B), are much more significant than for E~ . The dissimilarities with the exact 
results are also more pronounced, and especially in the former example. We do not find 
out any immediate explanation for this outcome. Yet, one should not be seriously worried 
by the results obtained in the case (A), as in this example there is only one proton in the 
final state. On the other hand, in the case (C) all three QRPA approaches for E'^ are 
quite similar before the collapse of the QRPA, and reproduce reasonable well the exact 
results. It should also be pointed out that only the SCQRPA is capable to account for 
the minimum of the exact E~^ . 

5.3 Double beta decay matrix elements 

The matrix element for the /3/?2i/ decay can be cast in the form: 



Ef+E^ 
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where |0) and |0) stand, respectively, for the ground states in the initial (A^, Z) and final 
{N — 2, Z + 2) nuclei, with energies Eq and Eq, and 

Ef = Ei- Eo; E+ = Ei - Eo, (5.9) 

are the corresponding excitation energies from the ground states to the virtual states |0;). 

As said before, the matrix elements (0|T_|0;) obey the sum rule (|5.4|) with given 
by ( ^.2|) . On the other hand, the matrix elements (0;|T_|0) fulfill the relation: 

S--S+ = N-Z~A, (5.10) 

where 

S^ = Y.m\Tm'- (5.11) 

I 

Within the QRPA calculations performed here there is only one intermediate state. 
Thus 

where, in the approximate calculations, (Oi„t|T_|0) and E~ are given, respectively, by eqs. 



J) and dJ) 



E+ = LU - Xp + Xrr, (5.13) 

and 



(o|r_|o„,) = ^2nDp4yupVn + XvpUn). (5.i4) 

Here the barred quantities correspond to the final nucleus. 

Before proceeding, let us note that we have found numerically that, within the exact 
calculation, only one intermediate state contributes significantly. This state is the isobaric 
analog state (IAS) of the state |0) in the isospin symmetry limit, i.e. when s = 1 and 



Fpn = 0. Thus, the eq. ( 5.12 ) can also be used in the exact calculations for all practical 



purposes. Moreover, we can write 

J Q- Q+ 

Therefore, in the 0(5) model the matrix element for the l3l32u decay only depends on 5*", 
E~, 5'+ and E~^. The results for the first two quantities have been discussed already, and 
the results for the remaining two are shown in figures 7 and 8. It should be specified 
that the QRPA equations for the final systems collapse at significantly smaller values of 
s, than for the initial ones. Namely, now Scrit = 1.17, 1.13 and 1.04, for the cases (A), (B) 
and (C), respectively. Nevertheless, from the comparison of figures 2 and 7 it can easily 
be seen that is quite similar to S'^. In the same way from figures 6 and 8 one realizes 
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that and are alike to each other, although the discrepancies between the exact 
and the approximates results are somewhat smaller for the former. Finally, in figure 9 are 
presented the results for the moments Ai2u- Note that the sign of A42u is not physically 
relevant. Only for esthetic reasons we will choose it to be positive before becoming zero 
and negative afterwards. From the comparison with the exact calculations one can infer 
that the all three approximate calculations account for the general behavior of A42u, in 
the sense that they all cross the axis s close to the point where the exact calculations 
do it. They, however overestimate its magnitude, for s < Scnt- From eq. ( |5.15|) and 
the inspection of figures 2, 5, 7 and 8 it is easy to discover that the main reason for the 
discrepancies with the exact results comes from the differences in E~^. Close to the critical 
value of s and beyond the SCQRPA is definitively "better" than RQRPA. Yet, in the case 
(C) all three approximations gives quite satisfactory results. 



6 Summarizing Discussion 

In this work we again took up the presently very debated issue of extending the quasiparti- 
cle RPA in the description of the charge-exchange transitions, to account more accurately 
for the physics around the transition point, which seems to be placed precisely in the 
region where the particle-particle force takes its physical value. An adequate tool to go 
beyond the QRPA approach seems to be its self-consistent extension, i.e. the SCQRPA. 
This approach properly includes ground state correlations and the Pauli principle, and 
has recently shown interesting results for the pairing problem as well as for other 



physical systems |2^, |29[. Because of the complexity of the problem it seems appropriate 
that this method be first tested on exactly solvable model cases before large scale realistic 
calculations are performed. 

It is considered |2^, ^ that the schematic 0(5) model is the appropriate testing 
ground for the double beta decay, and that it mimics rather accurately what is happening 
in realistic model calculations. But, one should be cautioned in doing such an extrapola- 
tion. First, the 0(5) model is only able to describe the F transitions that play a marginal 
role in the (3/3 decay, in comparison with the GT transitions. Second, the F processes are 
related with the isospin conservation that is a good symmetry in nuclear physics, while 



the GT processes are associated with the Wigner SU(4) symmetry ||3^, 0, which is badly 
broken by the mean field and only partially restored by the residual interaction. Third, 
in realistic calculations the QRPA collapses in the physical region for the particle-particle 
strength only in the GT case |P, |^ . 

To study specifically the GT transitions within a schematic model one has to resort 
to the more general S0(8) model |3^ where the T = and T = 1 pairings coexist. 
Recently it was shown that the RQRPA, when applied to this model, fails in re- 
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producing the beta and double beta decay properties beyond the phase transition point 
(SU(4) hmit). This failure was attributed to both: i) the need of changing the BCS basis, 
due to the transition to a T = pairing region, and ii) the violation of the Ikeda sum 
rule. One might hope that the SCQRPA would solve the problem. Yet, although SO (8) 
is only slightly more complicated than 0(5), we do not know how to build up explicitly 
the correlated ground state for the former, and, as was shown in the present work, this 
is the key point in solving the SCQRPA equation. The determination of the correlated 
ground state is also the main difficulty in applying the SCQRPA to real nuclei. Therefore, 
from the practical point of view, one of the most important tasks to solve the SCQRPA 
equations is the development of physically sound and numerically accurate approximate 
methods, which do not require the explicit construction of the correlated ground state. 
Schematic models such as the 0(5) and SO (8) are perfect testing grounds for this trial. 

In the present paper we have studied in detail the properties of the QRPA, RQRPA and 
SCQRPA equations in the 0(5) model. Compared with the previous work [^, we have 
made several improvements, namely: i) we have solved exactly the SCQRPA equations 
of the 0(5) model by evaluating everything exactly with the QRPA vacuum state ( ^.3.2| ) 
ii) we correctly took into account the coupling of the mean-field to the excitations. We 
also have examined carefully the RRPA limit of the model. The importance of all these 
features is clearly born out by the results. The energies, transition strengths, and double 
beta decay matrix elements are accounted for reasonably well. It clearly turns out that the 
SCQRPA approach performs better than either standard QRPA or renormalized QRPA 
for the description of these observables. Only in the case (C) the QRPA collapses in the 
physical region for the parameter s (s = 1). Thus, the samples (A) and (B) mimic better 
the F double beta decays and the example (C) the GT ones. 

Before the collapse of the QRPA all three approaches reproduce well the exact results. 
Near the transition point only the SCQRPA values are close to the exact ones. Beyond 
that point both the SCQRPA and RQRPA yield different results from the exact ones, 
but the former are always qualitatively superior. It can be suspected that this situation 
even prevails to some extent in the realistic situation. On the other hand, it is very well 
known that after the point where standard RPA breaks down, one has to change 
the single particle basis. In the present situation, where the proton-neutron correlations 
become "soft", this implies that after the phase transition point one has to consider the 
quasi-particles as mixed proton-neutron states. 

Elaborating a SCQRPA in this new basis should allow to go in a smooth way through 
the phase transition (as happens in other models ^2[)' expect that then the 

SCQRPA will reproduce the exact results not only before but also beyond the transition 
point. In the same line one then also should include proton-proton and neutron-neutron 
pair excitation terms in the QRPA operator using the generalized Bogoljubov technique. 
Unfortunately such a generalization of the theory is a nontrivial task even in the 0(5) 
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model and demands considerable formal developments. We intend to elaborate on this in 
the future. 

A novel feature that has emerged from this work is that in the RRPA limit not only 
the proton-neutron components of the interaction are renormalized, as was well known 
previously Ijl^, but also the quasi-particle energies. To get this effect it is essential to 
consider the coupling of the excitation to the mean-field, embodied in the generalized 
BCS equations . 

One further technical aspect that deserves mentioning here concerns the inclusion of 
so-called anomalous terms a^a in the excitation operator, as this was proposed in |^3|, |24 . 



Inclusion of such terms is very important to fulfill both the energy weighted and energy 
non weighted sum rules, as this was shown in for the case of homogeneous infinite 
nuclear matter. A straightforward calculation also shows that such terms are essential to 
fulfill the sum rule for F and GT transitions in realistic finite nuclei. Inclusion of these 
terms is also relevant for the correct treatment of symmetries and the appearance of the 



Golstone mode [^. But, contrary to the infinite matter case, the consistent inclusion of 
such terms in group theoretical models, like the one treated here, is a delicate problem 
and demands further studies. 

A particular feature of our model relates to the fact that the QRPA ground state wave 
function could be constructed explicitly and therefore the system of equations could be 
closed without further approximations. In the realistic situation it is certainly far from 
obvious to obtain the ground state wave function and approximations must be used to 
close the system. In this respect efficient methods have been developed in the past, like 
e.g. the number operator method [^]. On the other hand, it has been shown that the 
system of equations can be closed without knowing the wave function in coupling the 
SCRPA to the one body Green's function |]23[. This way was considered in eq. (16) of ref. 
2^ . In the future the elaboration of a more consistent single particle equation coupled 
to SCRPA will be a key point in closing the equations also for realistic situations. 

This work was supported in part by the DGICYT (Spain) under contract No. PB 
95/0123-C02-01. 
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Appendix A: Expectation values within the SCQRPA 

Except for the operators A^ppAl^,^, A^ppApp and AfnpJ^pn, the calculation of the expecta- 
tion values can be easily done using the completeness properties of the SCQRPA exchange 
mode and the properties of the basis states. For these three operators we had to use the 
generator state technique from ref. p9|| . 

Using the notation, 

{^^) = {o\^^pp\o) = {o\Kn\o), (a.i) 

and the ground state given by ( p.3.2|) we get 

n 

{M) = 2NoY.la^z^\2l\2l), (A.2) 

1=0 



(0|AA^,|0) = (OIACJO) = {0\K^Afpp\0) = 4Ar,E/2^2^2.^2/|2/), (A.3) 



n 

r 

1=0 



AM ^ 

mlpALlO) = -^^5:/(2fi-2/ + l)a,V'-i(2/|2/), (A.4) 



1=1 



mlpApplO) = (0|^LA„,|0) = ^^E^(2^ - 1)«,V'(2/|20, (A.5) 



1=1 



mlApnlO) = ^^{Q - (AT)), (A.6) 



and 



{0\ApnApn\0) = (OI^^^JO) = - m, (A.7) 

{0\AfnjMpn\0) = {0\AfpnAfnp\0) = 2{M) - {Q\AlpApp\Q) . (A.8) 
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Appendix B: RQRPA limit for the expectation values 

With (A/") defined in ( |A. 1| ) we obtain 



{o\Af;^\o) = m^jo) = {Mf + {Af) 



2n 



{0\AlA,,\0) = (Ol^LAnlO) = 2{0\AlA,n\0) 



and 



n 



2n 



2n 



The following relations are also useful: 



(2fi - Aff |0) = + ^{^n - l){{Af) - 2n), 



n 



{2n-i){{Af) -2n), 



{0\2{n-Afpp)'-AlApp\0) = 2in' + {Af){l-2n) 



(AO: 



(B.l) 

(B.2) 
(B.3) 

(B.4) 
(B.5) 

(B.6) 
(B.7) 



. (B. 
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Figure Captions 



Figure 1: Frequencies u for the {N,Z) systems for the cases (A), (B) and (C), as a 
function of the particle-particle coupling constant s. 

Figure 2: Transition strengths S~ for the {N,Z) systems for the cases (A), (B) and (C), 
as a function of the particle-particle coupling constant s. 

Figure 3: The ratio z = y /X (upper panel), expectation value of the number of quasipar- 
ticles in the ground state {J\f) = (0|A/^|0) = (0|A/'„|0) (middle panel), and the occupation 
probability for neutrons v"^ (lower panel). Within the QRPA, (A/") = 2^^. 

Figure 4: Excitation energies E~ for the case (B), as a function of the particle-particle 
coupling constant s. The approximate results were obtained with the prescriptions (|5.5|) 
(upper panel) and ( |^ ) (lower panel). 

Figure 5: Excitation energies E' for the cases (A), (B) and (C), as a function of the 
particle-particle coupling constant s. The approximate results were obtained from eq. 

Figure 6: Excitation energies for the cases (A), (B) and (C), as a function of the 
particle-particle coupling constant s. The approximate results were obtained from eq. 

&■ 

Figure 7: Transition strengths for the (A^ — 2, Z + 2) systems for the cases (A), (B) 
and (C), as a function of the particle-particle coupling constant s. 

Figure 8: Excitation energies E^ for the cases (A), (B) and (C), as a function of the 
particle-particle coupling constant s. The approximate results were obtained from eq. 

Figure 9: Matrix element Ai2u as a function of the particle-particle coupling constant s. 
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